function d=distance(i,j)
pathdir='mesh/';
path{1}=num2str(i);path{2}=num2str(j);
% [vertices1, faces1 ]= obj__read( [pathdir,'mesh',path1,'.obj'] );
% [vertices2, faces2 ]= obj__read( [pathdir,'mesh',path2,'.obj'] );
% P=vertices1';Q=vertices2';
% tic
parfor i=1:2
    vertices{i}=readobj_vertex([pathdir,'mesh',path{i},'.obj']);
end
% toc
% tic
[d, D] = HausdorffDist_par(vertices{1},vertices{2},1);
% toc

% tic
% pathdir='mesh/';
% path1=num2str(i);path2=num2str(j);
% % [vertices1, faces1 ]= obj__read( [pathdir,'mesh',path1,'.obj'] );
% % [vertices2, faces2 ]= obj__read( [pathdir,'mesh',path2,'.obj'] );
% % P=vertices1';Q=vertices2';
% [vertices1]=readobj_vertex([pathdir,'mesh',path1,'.obj']);
% [vertices2]=readobj_vertex([pathdir,'mesh',path2,'.obj']);
% P=vertices1;Q=vertices2;
% [d, D] = HausdorffDist_par(P,Q,1);
% toc